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Abstract 


In this paper, we describe a symbolic synthesis method which given 
an algebraic expression that specifies a bitstream function f, con- 
structs a (minimal) Mealy machine that realises f. The synthesis 
algorithm can be seen as an analogue of Brzozowski’s construction of 
a finite deterministic automaton from a regular expression. It is based 
on a coinductive characterisation of the operators of 2-adic arithmetic 
in terms of stream differential equations. 


1 Introduction 


A (binary) Mealy machine is a deterministic automaton which in each step 
reads an input bit, produces an output bit and moves to a next state. The 
induced mapping of input streams to output streams is a causal bitstream 
function, which we call the bitstream function realised by the Mealy ma- 
chine. In this paper, we describe a symbolic synthesis method which given 
an algebraic expression that specifies a bitstream function f, constructs a 
minimal Mealy machine that realises f. 

The inputs to our synthesis algorithm are called function expressions 
and they define bitstream functions in the algebra of 2-adic numbers. Here 
we use that the formal power series representation of a 2-adic integer can 
be seen as the bitstream of its coefficients. We describe the 2-adic algebra 
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below, but for now such a function expression can be thought of as a spec- 
ification of a function on the rational numbers. The interesting property of 
2-adic arithmetic is that it allows us to calculate with bitstream represen- 
tations of rational numbers in an easy manner similar to how one computes 
with integers. 

The synthesis algorithm described here can be seen as an analogue of 
Brzozowski’s construction in [3] of a finite deterministic automaton from a 
regular expression. In particular, we show that the set of function expres- 
sions carries the structure of a Mealy machine by giving an inductive defini- 
tion of derivative and output of function expressions. This Mealy machine 
of expressions is defined in such a way that the algebraic semantics coincides 
with the behavioural semantics of Mealy machines, which will be defined in 
terms of causal stream functions. Now given a function expression F speci- 
fying a function f, we obtain a realisation of f by the symbolic computation 
of the (sub)machine generated by F. (The generated submachine is in gen- 
eral not minimal, but we can ensure minimality by reducing expressions to 
normal form.) The language of 2-adic arithmetic allows the specification 
of functions that cannot be realised by any finite Mealy machine. But we 
shall identify a subclass of so-called rational function expressions for which 
a finite realisation exists. 

In the design of digital hardware, Mealy machines specify the behaviour 
of sequential circuits, and there exist algorithms which construct from a 
(finite) Mealy machine, a sequential circuit which exhibits the specified be- 
haviour. Combining these algorithms with our synthesis algorithm we thus 
obtain a complete construction from algebraic specification to sequential 
circuit. 

In summary, the main contributions of the present paper are: (1) The 
elementary but useful observation that the set of all causal stream functions 
constitutes a final Mealy machine, which forms the basis for a behavioural 
semantics of Mealy machines in terms of their minimisation. (Although 
minimisation of Mealy automata is well known [4], its characterisation here 
by means of finality is new.) (2) The insight that given a rational func- 
tion expression F, we can effectively construct by symbolic computation a 
(minimal) Mealy machine that realises the bitstream function specified by F. 

The basis for the present synthesis algorithm was given in [18]. These 
ideas were developed into a proper algorithm in [8] (for an implementa- 
tion, see [7]). Further results on complexity and size of realisations were 
included in the first author’s PhD thesis [6, Ch. 3]. This paper contains a 
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short, but improved presentation of the basic results in the abovementioned 
work. In particular, the presentation of the Mealy machine of expressions, 
the algebraic semantics of function expressions, and the proof that algebraic 
semantics coincides with behavioural semantics for function expressions are 
new with respect to [6, 8, 18]. Moreover, the synthesis algorithm and com- 
plexity analysis have been simplified by making use of particular properties 
of rational function expressions. Finally, we mention that the theory un- 
derlying the present work is essentially coalgebraic (cf. [15]), but we have 
deliberately chosen for a presentation which does not require any familiarity 
with coalgebra. 

Structure of the paper: In Section 2 we give the basic definitions re- 
garding Mealy machines. In Section 3 we describe the 2-adic algebra of 
bitstreams and define the function expressions which serve as our specifica- 
tion language. In Section 4 we show how function expressions can be turned 
into a Mealy machine, and in Section 5 we present the synthesis algorithm 
for rational function expressions and analyse its time complexity. Finally, 
we discuss related and future work in Section 6. 


2 Mealy Machines 


We give the basic definitions on Mealy machines, streams and causal stream 
functions. We introduce the notion of stream function derivative and show 
how it can be used to turn the set of causal stream functions into a final 
Mealy machine, thus providing a characterisation of minimal Mealy ma- 
chines. 

Let A and B be arbitrary sets. A Mealy machine (S,m) with inputs 
in A and outputs in B consists of a set of states S and a transition function 


m: S > (Bx S)4 


This function maps a state so € S to a function m(so): A > (Bx S), which 
produces for every input a € A a pair (b, 51), consisting of the output b and 
the next state s;. We call a Mealy machine binary if inputs and outputs 
are taken from the set 2 = {0,1}. We will write 


alb 


s t iff m(s)(a) = (b,t). 


For an arbitrary set A, we denote by A” the set of streams over A, i.e., 
AY = {a | a: N > A}. An element a € A” will be denoted by a = 
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(a(0), a(1), a(2),...). The head and tail maps on streams are denoted by 
hd and tl, respectively, that is, for a € AY’, hd(a) = a(0) and tl(a) = 
(a(1), a(2), a(3),...). Later, in the context of stream differential equations, 
we will also refer to hd(q@) and tl(q) as the initial value and stream derivative 
of a, respectively, and write a’ instead of tl(a@). Moreover, for a € AY and 
n €N, we write al, for the finite prefix (a(0),...,a(n)). 

We define the (input-output) behaviour of a state so in S = (S,m) 
as the stream function behs(so): AY — BY which maps an input stream 
(ao, @1,a2,...) € A” to the output stream (bo, b1, b2,...) € BY given by the 
unique sequence of transitions 


ao|bo ay|by ap|be ap41|be41 
> S| Shs, cae > Sk4y + 
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If S is clear from the context, we will often leave out the subscript and simply 
write beh(s). We say that a state s in (S,m) realises a stream function f if 


beh(s) = f. 


Example 1 The figure below shows an example of a binary Mealy machine 
which starting in state sg counts the number of ones in the input modulo 
2. Formally, so realises the function f defined on input stream a € 2” by 
f(a)(n) = yep a(4) mod 2 for alln EN. 


0|0 o|1 1/0 0/0 
aS 
S: ay ay fo. 
ie Wee 
11 41 


Note that, for alla € 2” andn €N, f(a)(n) is determined by a(0),..., a(n). 


We call a stream function f: A” > BY” causal if the n-th element of 
the output stream f(a) depends only on the first n elements of the input 
stream. More formally, f is causal if for all n > 0 and all a, @ € A®: 


aln= Bln = > — fla)(n) = f(8)(m). 


It is straightforward to prove that for any Mealy machine (S,m) and any 
s € S, beh(s): AY > BY is causal. Let I denote the set of all causal stream 
functions, i.e., 


T={f|f: A* — BY and f is causal }. 
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Next, we will define a function y: Tl — (B x T)4 such that (I,y) is a 
Mealy machine. For a € A” and a € A, we denote by a: a the stream 
(a, a(0), a(1),a(2),...), and given f: A” > BY”, we write f(a:—) for the 
stream function that maps a to f(a: a). Now let f € T anda € 2. We 
define 
fla] := hdof(a:-) € AX >B (1) 
fa i= tlof(a:—-) € AY” BY 


Since f is causal, it follows that also fa is causal (hence fa € I’), and that 
fla] is constant, hence f[a] can be considered an element of B. We call f|a] 
the initial output of f on input a, and fa is the stream function derivative 
of f on input a. We define the transition function 


qr TS(BxT)4 by y(f)(a) = (fla), fa) for all f ET,ae A. 
This gives us an (infinite) Mealy machine [ = (I, y) with transitions 


al fla] 
{Ja 
Next we characterise (I, y) using the following notion. A homomorphism 
of Mealy machines from (S,m) to (S’,m’) is a function h: S — S’ that 
preserves transitions: if m(s)(a) = (b,t) then m’(h(s))(a) = (b, A(t)); in 
other words, 


alb 


s— >t => hs) ge 


——> hit) 
Theorem 2 For any Mealy machine S = (S,m), the map behs: S >T is 


the unique homomorphism from S to T. In other words, T is a final Mealy 
machine. 


Proof: Let (S,m) be an arbitrary Mealy machine. We denote the output 
and next state functions defined by m at state s € S by og, and dg, re- 
spectively, that is, m(s)(a) = (0s(a),d5(a)). To see that beh: S >Tisa 
homomorphism, let s € S,a € A anda € A” be arbitrary. We have: 


beh(s)[a] = hd(beh(s)(a:a)) = 0;(a), 
and by letting so = ds(a) and s;41 = ds,(a(z)) for all 7 > 0, we have 


(beh(s)a)(a) = tl(beh(s)(a:a)) 

= t(0s(@), Os9(@(0)), 0s, (a(1)), -- +) 
(0s9(@(0)), 0s, (a(1)), -- +) 
beh(d;(a))(q). 
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Hence beh(s)q4 = beh(d;(a)). The proof that beh: S > T is the unique 
homomorphism from S to T is left to the reader. 


The universality of ([,y) can be expressed in yet another way. We 
need the following notions. For a state s € S of a Mealy machine (5, m), let 


(s)) CS 


denote the smallest subset that contains s and is closed under transitions 
(for any inputs). Clearly, ((s)) is also a submachine of (S,m), by taking 
as its transition function the restriction of m to the set ((s)). We call ((s)) 
the submachine generated by s in (S,m). We call ((s)) a realisation of f, if 
behs(s) = f, i.e., s realises f (in S). We say that a Mealy machine (S,m) is 
minimal if beh: (Sm) — (I, y) is injective. From the finality of T, we get: 


Corollary 1 For all f €T, (f)) is a minimal realisation of f. 


Proof: Follows from the fact that the inclusion map ((f)) — T is an injective 
homomorphism of Mealy machines. 


The final Mealy machine T thus consists of all the behaviours that can 
be realised by some Mealy machine, and we therefore refer to beh(s) as the 
behavioural semantics of a state s. 

For a causal stream function f: A” — BY”, we use the following nota- 
tion for repeated stream function derivatives: for w € A* and a € A, we 
define f- = f (where ¢ is the empty word) and fwa = (fw)a. Hence the state 
set of ((f)) equals the set {fy | w € A*} of all stream function derivatives of 


f. 


Example 3 We illustrate the notions and results above with the binary 
counter from Example 1. Let f = beh(so), g = beh(si) and h = beh(s2). It 
can easily be seen that f = h and g(a)(k) =1— f(a)(k) for all a € 2” and 
k € N. Moreover, we have the following initial outputs and stream function 
derivatives: 


f[0] = 0, ffl] = 1, g|0] = 1, g|1] = 9, fo=f, fi =49; Jo = 9; n=f 
and thus beh maps S = ((so)) to its minimisation ((f)) C T, on the right: 


0|0 O]1 10 0|0 0|0 11 O|1 
1] 1|0 
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In algebra, the behaviour of Mealy machines is typically described in 
terms of functions of type At + B. Although this set is isomorphic to I, 
we prefer to work with the latter because of the rich algebraic structure on 
streams. The notion of stream function derivative already occurs in [13], 
and is called state there. It is also a variation on the classical notion of 
derivative (or inverse) of functions from A* to B* (cf. [4]). Also Theorem 2 
and Corollary 1 are essentially reformulations of classical results. 

The main contribution of the present paper consists of the observation 
that ((f)) can often be constructed by a symbolic computation of stream 
function derivatives starting from an algebraic specification of f. In Sec- 
tion 5, we shall describe such a symbolic algorithm for bitstream functions 
specified in 2-adic arithmetic. 


3 The 2-Adic Bitstream Algebra 


We will specify bitstream functions in the algebra of 2-adic integers (cf. [5]). 
A 2-adic integer is usually written as a (formal) power series of the form 
> p29 42’ where a; € 2 for all i € {0,1,2,...}. We identify such a power 
series with the bitstream (ag, a1, @2,...), so for us the set of 2-adic integers 
is simply the set 2” of bitstreams. 

There is a (strict) inclusion of the set of rational numbers with odd 
denominator 


Qoaa = {p/a | p,4 € Z,q odd } 


into the 2-adic integers by taking infinitary base 2 expansions (see e.g. [10]). 
For a positive integer n, this is just the binary representation of n (least 
significant bit on the left) padded with a tail of zeros; for instance, 


Bin(2) = (0,1,0,0,0,...) = 010" 
Bin(5) = (1,0,1,0,0,...) = 1010” 


Binary representations of negative integers end with an infinite sequence of 
ones and rational numbers (with odd denominator) have binary representa- 
tions that are eventually periodical; for instance, 


Bin(—1) = (1,1,1,1,...)=1% 
Bin(—5) = (1,1,0,1,1,1,...) = 1101” 
Bin(1/5) = (1,0,1,1,0,0,1,1,0,0,1,1,0,...) = 1(0110)”. 


(Rationals with even denominator require formal power series representa- 
tions )>°°, a;2' where k < 0, and are not treated here.) 
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Below we shall define the binary representation Bin(q) of any rational 
gd € Qoaa by means of a stream differential equation. Such equations specify 
streams o, in analogy with traditional differential calculus, in terms of their 
initial value o(0) and derivative o’ (which are defined as hd(c) and tl(c)). 
For instance, the differential equation 


=o -o0y=1 


clearly defines the stream (1,1,1,...). We refer to [16] for an overview on 
stream differential equations. 
Now let odd(n/2m +1) = n mod 2 for n,m € Z. We define the 
inclusion map 
Bin: Qoda — 2” 


by the following system of stream differential equations (one for each q € 


Qoaa): 
Bin(q)(0) = odd(q), Bin(q)’ = Bin((q — odd(q))/2) (2) 


The reader is invited to compute the examples of binary expansions given 
above, using this definition of Bin. 

We recall that the set of rational numbers (with odd denominator) is 
an integral domain: a commutative ring in which multiplication has no zero- 
divisors, i.e., for all p and q, if p x gq = 0 then p = 0 or q = 0. Next we shall 
introduce operations of addition and multiplication (as well as minus and 
inverse) on the set of 2-adic integers such that they reflect the operations 
on the rationals. More formally, we shall turn the set of 2-adic integers into 
an integral domain 


Aovadic = (2°, +, 7; Xs [0], [1]) 


such that the inclusion map Bin : Qogg > 2” is a homomorphism of integral 
domains. 

In the literature, the operators on the 2-adic integers are often defined 
by explicitly specifying how they work on binary representations, see Fig- 
ure 1 for examples of addition and multiplication. Their definitions are 
sometimes also given in terms of their representation as power series, typi- 
cally in the form of some recurrence relation on their coefficients. Here we 
shall give, instead, a definition by means of stream differential equations, 
which we saw already examples of above, and which can be seen as a gen- 
eralisation of definitions by recurrence, cf. [16]. Our choice to use stream 
differential equations is not just a matter of taste: they form the basis for 
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addition: Lilies (carry) 
101000--- = 9) 

+011111--- = —2 

110000--- = 3 


multiplication: 11000---x O1111--- = 3x (-2) 
O00000--- 
11000.-.-- 
11000--- 
11000--- 
11000--- 


+ 


O10111--- = -6 
Figure 1: Examples of 2-adic addition and multiplication. 


the definition of the Mealy machine of stream function expressions, to be 
introduced in Section 4. 
First we define the constants zero and one by 


[0] = 00,0705 2) (2) E0005 -.2) 


Below we shall use the following Boolean operators on 2 = {0,1}, which 
are defined, for all a,b € 2, as usual: a A b = min{a,b} and a@b = 1 iff 
a=0,b=1 ora=1, b=0 (exclusive or). 


Definition 1 We define the operators of addition, minus, multiplication 
and inverse of 2-adic integers by the following system of stream differential 
equations. For a, B € 2”, 


| derivative: initial value: 
(a+ BY = (a’ +") +[a(0) A BO) | (a+ 6)(0) = a(0) 6 B(0) 
( —(a! + [a(0)}) (—a)(0) = a(0) 
(ax B)! = (a’x B) + ([a(0)] x 8’) | (ax B)(0) = a(0) A B(0) 
( —(a! x (1/a) ) (1/a)(0) = 1 


There is a side condition for the definition of inverse: 1/a is defined only 
for a with a(0) = 1. 


| 
o 
| 


be 

™_ 
Q 

Rac 
| 
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We briefly explain the intuition behind the equations above. (For the 
fact that the operators above are uniquely defined by their defining stream 
differential equations, we refer to [16].) The equation for sum shows that 
a carry term must be added in case the two initial values are both 1. The 
equation for minus is obtained from the requirement that (—a) +a = [0] for 
alla € 2”. By taking initial value and derivative on both sides and using the 
equation for +, we find that (—a)(0)®a(0) = 0, hence (—a)(0) = a(0), and 
that (—a)’ + a’ + [a(0)] = [0], hence (—a)/ = —(a’ + [a(0)]). The equation 
for the product states that for all a, 8 € 2”, a x 6 can be calculated using 
the base 2 version of shift-add-multiplication known from the multiplication 
in decimal notation. Finally, the multiplicative inverse of a € 2” is defined 
only if a(0) = 1, since there is no bit a € 2 such that aA a(0) = 1. If 1/a 
is defined then it satisfies (1/a) x a = [1]. The equation for 1/a can be 
derived in a similar way as for —a. 

The operations on the 2-adic numbers have been devised in such a way 
that Bin: Qoaq 4 Agadic is a homomorphism of integral domains: one can 
easily show that, for all p,q € Qoaa; 


Bin(p x q) = Bin(p) x Bin(q) 


(with on the left multiplication of rationals and on the right multiplication 
of bitstreams), and similarly for the other operators. 

When calculating with the 2-adic operations it is convenient to have a 
constant denoting the bitstream Bin(2). We define 


X := Bin(2) = (0,1,0,0,0,...), 


The constant X can be used to express some identities on bitstreams that 
will be useful later: For all a € 2”, 


[a(0)] + (X x a’) 
ata = Xxa 


. (3) 
4 Mealy Machine of Expressions 


We will specify bitstream functions in the language of Aggie over a single 
variable s. Formally, the set of function expressions FExpr is generated by 
the following grammar: 


FG s:= s|O0|1|X|—-F|F+G|FxG|1/(1+(X x F)) 
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We use the symbols —,+,x,/ to denote the operations on bitstreams as 
well as the corresponding syntax constructors. The typing should always 
be clear from the context. We will use standard notational conventions: we 
write F” for the n-fold product of F with itself, (in particular, Fo = 1) and 
F/(1+(X x G)) or a) instead of F x (1/(1+ (X x G))). The algebraic 
semantics of function expressions is given by the expected interpretation in 
Aoadic. The reason we only allow division by terms of the form 1 + (X x F) 
is to ensure that the inverse operation is always defined when evaluating 
function expressions in Aagdic. 


Definition 2 We define the algebraic 2-adic semantics [F]: 2” > 2” of a 
function expression F € FExpr by the following inductive clauses. Let o be 
a bitstream, 


=o ©: [-Fl(o 


( 
0] [F + G]( 
( 
( 


II 
an 
x 


1] [F x G] 
ss [1/(1 + (X x FI (o 


) 
) 
g) 
) 

We say that a function expression F specifies the bitstream function 
[F], and two function expressions F and G are equivalent (notation: F = G) 
if [F] = [G]. 

One easily shows (by structural induction on function expressions) that 
for all F € FExpr, [F] is a causal bitstream function. In other words, 
[—] is a map from FExpr to T. We will now define a transition function 
€: FExpr — (2 x FExpr)? such that (FExpr,€) is a binary Mealy machine 
and the algebraic semantics [—]: FExpr — [ is a Mealy homomorphism. 
That is, for each F € FExpr and a © 2, we want to define Fla] € 2 and 
F, € FExpr such that Fla] = [F][a] and [Fa] = [F]a. In order to do so, first 
recall how we defined 7(f)(a) in terms of f(a:—), hd and tl (cf. equation 
(1)) for f €T and a € 2. We will “mimic” y in the syntax. 

First, we need a syntactic version of the map f(a:—) for f € T and 
a € 2. More precisely, given F € FExpr and a € 2, we want to find a 
function expression which specifies [F](a:—). Note that for all o € 2%, 
[0 + (X x s)](o) =0:0 and [1+ (X x s)](o) = Lo. 


Definition 3 We define 1: 2 + FExpr by 1(0) = 0 and u(1) = 1, and for 


F € FExpr anda € 2, we let a:s := u(a) + (X x s). We define F(a:s) to be 
the function expression obtained by uniformly substituting a:s for s in F. 
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We now show that F(a:s) indeed specifies [F](a:—). 


Lemma 1 For all F € FExpr, alla € 2 and allo € 2”: 
[F(a:s)](o) = [F](a:¢). 


Proof: The lemma is proved by a straightforward induction on the structure 
of F. We only show a few example cases: 


[s(a:s)](c) 


(cf. eqn. (3)) 


[(F + G)(a:s))](c) 


u(a) + (X x s)](o) 
La) (a) + (EXC) = [s](o)) 


The transition function € on function expressions is defined inductively 
over the syntactic structure. In order to motivate the definition of €, con- 
sider, for example, a product expression F x G, and suppose that [F] = f 
and [G] = g. We then want [(F x G)a](o) = tl(f(a:o) x g(a:o)) for all 
o € 2”. By the stream differential equation for product, this means that 


[(F x G)a]() 


= (fla 


0) 
= (fala) x g 
yx 


= ([F]a(c 


) x gla:o)) + ([hd(f(a:e))] x (g(a:))) 
(a:o)) + ([flal] x ga(o) 
[s](a:o)) + CFI lal] x [¢]a(o). 


Comparing this equation with the definition of (F x G)q below it is easy to 
see that an induction argument and Lemma 1 will show that [(F x G)a] = 


[F x Gla. 
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Definition 4 Let €: FExpr > (2 x FExpr)? be the map given by €(F)(a) = 


(Fla],Fa) where Fla] and F, are defined for all F € FExpr anda € 2 by: 
syntactic initial output 
sa] = «| Fla = Flal 
Ola] = 0]| (F+G){a] = Fla] 6 Gla] 
1a] = 1] (Fx G)[{a] = Fla] \Gla] 
Xa} = 0] 4/4+(xF)))lal = 1 
syntactic stream function derivative 
Sa = 8s|| (-Fha = —(F, + (Fla])) 
0. = O| (F+G)q = (Fa+Gq) +4(F[a] A Gal) 
te, = OU SG) = (F, Xx G(a:s)) + (v(Fla]) x Ga) 
a = 1) G/4+(QXF)))a = —(Fla:s))/( + (Xx F(a:s))) 


We call (FExpr,€) the Mealy machine of expressions, and we refer to Fla] 
and F, as the syntactic initial output, respectively the syntactic stream func- 
tion derivative, of F on input a. 


We now show that € indeed ensures that the algebraic semantics con- 
cides with the behavioural semantics of function expressions. 


Proposition 1 The map [—]: FExpr — IT is a homomorphism of Mealy 
machines from (FExpr, €) to ([,7). 


Proof: We must show that for all F € FExpr and all a € 2: 


[F][e] =Fla] and [Fa] = [Fla. 


The proof is by induction on the structure of F. We only show the case 
for s and product; the others can be shown along the same lines (the full 
proof details are found in the Appendix). In the identities below, IH refers 
to the induction hypothesis. Note also that for any b € 2 and ao € 2”: 
[e(6)](o) = [b]. Let a € 2 and a € 2”. 


[s]lz] = ([s](a:c))(0) = (a:¢)(0)=a = sla]. 
[sla() = ([s](a:c))!) = (aio)! =o = [su](c). 
[Fx G][a] = ([F x G](a:c))(0) = ([F](a:¢) x [G](a:))(0) 
= [F]la) A [¢] [a] 
2 Fla] A Gla] = (F x G) [a] 


[Fx G]a(o) = 


Since (I,y) is a final Mealy machine, the map [—] coincides with the 
unique homomorphism beh: (FExpr,€) — (I, y). In other words, we have 
shown that for all function expressions F, [F] = beh(F). 


5 Synthesis 


A consequence of Proposition 1 is that the generated submachine ((F)) is 
a realisation of [F]. Conceptually, we can construct ((F)) by computing 
the transition closure of {F} in (FExpr,&). However, in general, (F)) is 
neither finite nor minimal. In order to obtain a minimal realisation of |F] we 
compute ((F)) modulo equivalence. In practice, this is achieved by reducing 
function expressions to normal form, which we briefly describe now. 


5.1 Normal forms 


We call a function expression F integral if it does not contain the inverse 
operation, and F is closed if it does not contain s. The polynomial normal 
form pnf(F) of an integral expression F is an analogue of the distributed 
normal form of polynomials (in the variable s). It can be computed in the 
expected manner by applying identities of commutative rings: (i) distribute 
x over +, (ii) reduce using identities for 0 and 1, (iii) collect terms on powers 
of s, and finally (iv) reduce sums of closed expressions using the identity 
X” +x" = x"+!. for all n > 0. This last identity holds since for all a € 2”, 
ata=X xa, cf. equation (3). For example, pnf((1+X) x (1+s)4+1)= 
X? + (1 +X) x s, and this polynomial normal form is “computed” in the 
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following series of identities: 


((14+X)x(1+s))+1 (1+X)x1)+((1+X)xs)+1 

(1x 1)+(Xx 1))4+ ran (Xxs))+1 
1+X)+(s oe 
1+1+X)+(1+X)x 


X?+(1+X)xs 


~~ ~~ 


Note that in the last step we used that, 
fac yy See See ary 


For any integral expression F, it should be clear that F = pnf(F) and that 
pnf(F) is unique: pnf(F) = pnf(G) if and only if F = G. Given arbitrary 
function expressions F and G, we can therefore decide whether F = G by first 
rewriting (using the identities of integral domains) F and G into fractions 
P/Q and R/S, respectively, where P,Q,R,S are integral, and then checking 
whether pnf(P x S) = pnf(R x Q). These normal forms are treated in more 
detail in [6]. In fact, in our synthesis algorithm we only need to compute 
normal forms of closed, integral expressions, as we will see in Section 5.3. 


Example 4 We illustrate by computing (a representation of) (([F])) for the 
function expression F = (1+ X) xs. For the transition on input 1 we find 
that: 


((1+X)xs)[l] = (1{Jex{l]) As{l]=(e@0)A1=1. 
(14%) xs): = ((1+X), x s(1:s)) + (1 + XL) x 81) 
= (((4r + 4X1) +e(4[1] AX[I)) x (1 + (Xx 8))) 


+((1[1] @ X[1]) x s1) 
= ((0+1)+(e(1A0) x (14+ (Xx s)))) + ((1 60) x s) 
= (((0+1)4+0)x (1+(Xxs))+(1 xs) 
= 14+((14+X)xs)=14+F 


Computing further derivatives and initial output, we find the following 
minimal realisation of [F]: 


0/0 aap 10 11 
aa er 
ii fi fu 
oll 0/0 


where f = [F], f: = [1+ F] and fi: = [X+F]. 
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The normal forms essentially allow us to compute submachines in the 
final Mealy machine. Still, for arbitrary F € FExpr, the least fixed point 
construction of (([F])) is not guaranteed to terminate as it is easy to specify 
functions that have no finite realisation. For example, if we take F=s xs 
and we compute the derivatives [F]o, [Floo, [Fllooo,... we get the sequence 
[X x F], [X? x F], [X° x F],... which are all distinct. A similar argument 
shows that if F = 1/(1+(Xx-s)), then [F] has infinitely many distinct 
stream function derivatives. 


5.2 Rational functions 


We now define a class of expressions that specify functions with finite re- 
alisations. A rational function expression is a function expression of the 


form 
_D+(Cxs) 


~~ 1+(Xx&£) 
where D,C,E € FExpr are closed, integral function expressions. In other 
words, D,C and E specify constant bitstream functions whose value is of the 
form Bin(a) for some integer x € Z. An example of a rational function 
expression is 
(—X) + ((1+X) xs) 
1+ (X x (X+X)) 


which specifies the bitstream function 


_ Bin(—2) + (Bin(1+ 2) xo) _ Bin(—2) + (Bin(3) x o)) 
ie) —BadsaxO+D) Bin(9) 


(4) 


A function f: 2” — 2” is called rational if there is a rational function 
expression F such that [F] = f. Hence the function f(o) = Bin(3) x o 
specified in Example 4 is also rational by taking D=E=0 andC=1+4X. 
The numeric interpretation of closed, integral expressions will be convenient 
below. 


Lemma 2 A function f: 2” > 2” is rational iff there are d,m,n € Z such 
that n is odd, and for allo € 2” 


= Bin(d) + (Bin(m) x o) 


fo) Bin(n) 


Proof: This follows essentially from the fact that Bin: Qoag > Aaadic is a 
homomorphism of integral domains. 
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To simplify notation, we will leave out the Bin-part, and just write 
f(o) = (d+ (mx o))/n whenever f(c) = (Bin(d) + (Bin(m) x @)))/Bin(n). 
Similarly, for x € Z, we write x(0) and 2’ instead of Bin(x)(0) and Bin(z)’, 
respectively. The following technical lemmas will be used to prove that 
rational functions have finite realisations. Their proofs can be found in the 
Appendix. The first of these lemmas uses the numeric interpretation of 
rational functions to characterise the immediate derivatives. 


Lemma 3 Let f be a rational bitstream function of the form: 


_ d+(mxo) 
a n 


f(o) 


for integers d,m andn with n odd. Fora € 2, the stream function derivative 
fa is given by: 
6(a)+(mxoa 
(fa)(o) = Atm xe) (5) 


where (in the numeric interpretation) 


5 d if d even 
6(0) = 3? | (1) = 
3(d—n) ifd odd 
Hence Lemma 3 already tells us that the derivatives of a rational func- 


tion are again rational. The next lemma uses the numeric interpretation to 
give a bound on the range of 6-values that can occur in the stream function 


(d+ my) if d+m even 
(d+m-—n) ifd+m odd 


NIP NI 


derivatives of a rational function. 


Lemma 4 Let f(o) = dt(mxe) be a rational function with n > 0 odd. For 
all w € 2*, the stream function derivative fy is of the form 


d(w) + (m x o) 


(fw)(o) = (6) 
where 6(w) is an integer such that 


min{d,-n+1,-n+m+1} < d(w) < max{d,m-—1,0}. 


The crucial properties that make it possible to perform synthesis from 
rational function specifications are stated in the following proposition. 
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Proposition 2 For all rational bitstream functions f: 2” > 2%, 


1. all stream function derivatives of f are rational, 
2. (f)) is finite. 


Proof: Let f(a) = dt(mxo) be a rational function. Item 1 of the proposition 


follows from Lemmas 2 and 3. To see that item 2 holds, first note that we can 
always assume that n > 0 since f is equal to the function g(a) = —dtemxe) 


The number of states in (f)) equals the number of distinct 6(w)-values in 
the derivatives of f. By Lemma 4 this number is finite. 


5.3 Algorithm 


We construct (([F])) for a rational function expression F by computing states 
and transitions in the final Mealy machine. Stream function derivatives of 
[F] are denoted by rational function expressions, and in order to efficiently 
determine when two expressions denote the same derivative, we normalise 
the closed integral subexpressions of F at the beginning of the computation. 
More precisely, if F = (D+ (C x s))/(41 + (X x E)) is a rational function 
expression, then we define the reduced form of F as 


red(F) = (pnf(D) + (pnf(C) x s))/(1 + (X x pnf(E))), 


and we call F reduced, if F = red(F). Starting from a reduced expression, 
all derivatives computed during the synthesis algorithm are represented in 
a similar, reduced form. This is achieved by using the function next. 


Definition 5 Given a rational function expression 


D+(Cxs) 
fo eee 
1+ (Xx E) 
we define 
nextp(F,0) = Do + —(D[0] x E), and 
nextp(F,1) = ((D1 +Ci) +e(D[1] A C[1])) + —((D[1] @ Cf1]) x B), 


and fora € 2 we define 


pnf (nextp(F, a)) + (C x s) 


next(F,a) = T+ (kx 5) 
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It should be clear that starting from a reduced function expression 
F, any two expressions computed with the next-function are equivalent if 
and only if they are (syntactically) equal. The next lemma shows that the 
next-function can be identified with the syntactic derivative function. 


Lemma 5 For all rational function expressions F and alla € 2, 
next(F,a) = Fg. 


Proof: The lemma follows easily by writing out the details of the definition 
of F,, and applying integral domain identities. See also the proof of Lemma 3 
(in the Appendix). 


The synthesis algorithm is shown in Figure 2 (and we explain it in some 
more detail below). We denote the empty list by [], the concatenation of 
two lists « and y by x.y, and we let remdup be a function that removes 
duplicates from a list. If S$ is a list of states (i.e. reduced rational function 
expressions) and T is a list of transitions, then transitions(S) is the list 
of transitions with source state in S, and targets(T) is the list of target 
states of transitions in T. Using list comprehension, we can write this more 
precisely as: 


0|¢(0] 161] 


transitions(S) := [| G—=+> next(G,0),G ——> next(G,1))|GeS], 
targets (T) = (c¢|oesHeT]; 

1. NewStates := [red(F)]; States := []; Trans := []; 

2. do until NewStates = { 

3 NewTrans :=_ transitions(NewStates); 

4 Trans := Trans.NewTrans; 

5. States := States.NewStates; 

6 NewStates := remdup(targets(NewTrans)) \ States; 

7. 3} 

8. return Trans; 


Figure 2: Synthesis algorithm 
The synthesis algorithm is a standard least fixed point algorithm. We 


represent a (partially constructed) Mealy machine as a list Trans of tran- 
sitions. A list States is also used to keep track of which states have been 
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processed, and a list NewStates contains the states that were newly found 
in the previous iteration. Initially, NewStates only contains the initial spec- 
ification in reduced form. In each iteration, we compute the list NewTrans 
of transitions starting from NewStates, and add these to Trans. The new 
states for the next iteration are the target states of NewTrans except for 
the ones that have already been processed. At the end of iteration number 
k, all states and transitions up to depth k have been created. Hence when 
no new states are found, Trans represents (([F])). 

We can now finally state our synthesis result for rational function spec- 
ifications. 


Theorem 5 For any rational function expression F, we can effectively con- 
struct a finite, minimal realisation of [F] using the algorithm in Figure 2. 


Proof: By Corollary 1 and Proposition 2, (([F])) is a finite, minimal re- 
alisation of [F]. The algorithm in Figure 2 constructs a Mealy machine 
(S,m) whose states are (rational) function expressions. Due to Lemma 5, 
and the fact that all rational function expressions in S are reduced and have 
the same “denominator” expression, two function expressions $1, 52 € S' are 
equivalent if and only if they are (syntactically) identical. Hence (S,m) is 
minimal. 


We illustrate our synthesis algorithm with a couple of examples. 


Example 6 We construct the minimal realisation of the function f speci- 
fied by F = X? x s, i.e., in numeric notation, f(a) = 8 xo for all o € 2%. For 
the sake of readability, we will write derivative expressions in their numeric 
interpretation. For example, the derivative (X? + (X° x s)) will be denoted 
4+ 80. The diagram in Figure 3 shows the Mealy machine obtained from 
F using our algorithm. At the beginning of the computation, the only state 
is 8a. In the first iteration, the immediate derivatives of 80 are computed. 
These are 80 and 4+ 80, so 4+ 8a is a new state. In the second iteration, 
we compute the immediate derivatives of 4+ 80 and find 2+ 80 and 6+ 80, 
both of which are new. In the third iteration, we find the new states 1+ 80, 
5+80, 3+80 and 7+8co. In the fourth iteration, we find that all derivatives 
of the states 1+ 80, 5+80, 3+80 and 7+8o have already been visited, hence 
there are no new states after this round, and the algorithm terminates. 
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Figure 3: Mealy machine constructed from F = X° x s. 


Example 7 As yet another example, in Figure 4 we give a minimal realisa- 
tion of the rational function from equation (4), i.e., f(7) = ((—2)+(3xa))/9. 
For a compact presentation, the states are labelled only by the 6-value. For 
example, the state which realises the function (1+ (3 x 7))/9 is labelled just 
with 1. 


Figure 4: Mealy machine realising f(7) = ((—2) + (3 x a))/9. 
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Knowing that rational functions have finite realisations, it is natural 
to ask whether the converse holds, that is, whether all causal bitstream 
functions that are realised by a finite Mealy machine can be specified by a 
rational function expression. This question is answered in the negative by 
the following example. 


Example 8 Consider the Mealy machine depicted in the following diagram: 


oll 0|0,1|1 


onme 


We will show that there is no rational function expression F such that [F] = 
beh(q). To this end, we first observe that the behaviour of the state s is 
simply the identity function on bitstreams, hence beh(s) = [s]. Suppose, for 
the sake of arriving at a contradiction, that F is a function expression such 
that [F] = beh(q), and F is of the form F = (B+ (C x s))/(1 + (X x E)) where 
B,C,E © FExpr are closed, integral function expressions. By Proposition 1, 
[Fo] = [Flo = beh(q)o = beh(s) = [s], hence Fp = s. On the other hand, 
Lemma 3 tells us that Fo = (D+(C x s))/(1+(X x E)) for some closed, 
integral function expression D, and hence (D+ (C x s))/(1+(X XE)) =s. 
Consequently, D = 0, C = 1+(X x E), which implies that C[0] = C[1] = 1. 
From Definition 4 of the Mealy machine of expressions it follows that for all 
ae 2: 


Fla] = (B+(Cxs))[a] = Bla] 6 (Cla]Asla]) = Bla] @ (Aa) 


and hence, 
F/O] = B[O] and Fil] = 1-B{1] 


But since B is a closed function expression, B/0] = B[1], and hence by the 
above F[0] # F{1]. But beh(q)[0] = beh(q)[1] = 1 which contradicts the 
assumption that [F] = beh(q), and we conclude that no such F exists. 


Given the numeric nature of rational function expressions and their 
high level of abstraction, we do not find it surprising that they cannot 
specify all finite state Mealy behaviours. 
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5.4 Complexity 


In this section, we analyse the time complexity of our synthesis algorithm. 
Given a rational function expression F, we denote the Mealy machine con- 
structed by the synthesis algorithm by (S,mp) . We also need the fol- 
lowing definitions. The length Jen(E) of an integral function expression 
E is the number of symbol occurrences in E. In particular, len(s) = 1, 
len(E+G) = 1+ len(E) + len(G), and len(X*) = 2k — 1. If E is a closed, 
integral function expression, then val(E) is the integer « € Z such that 
[E] = Bin(x). For example, val(X? +1) = 5. For an integer x € Z, we 
denote by C, the unique closed, integral function expression in polynomial 
normal form such that val(C,) = x. 

We consider the computation of red(F) from F in line 1 of the algo- 
rithm as a preprocessing step. Consequently, we ignore its time cost in the 
overall analysis, and in the rest of this section, we assume that the initial 
specification is already reduced and of the form 


= Cg + com x s) 
Ch 


F (7) 
for integers d,m,n € Z with n > 0 odd. In particular, [F] is the function 
fo) = eEUENES. Moreover, we let K := |d| + |m| + |nl. 

As a first step, we give a high-level description of the complexity in 
terms of the subcomponents of the algorithm. 


Proposition 3 The synthesis algorithm which constructs from initial spec- 
ification F the Mealy machine (Sp,mg) runs in time O(RM + EM?) where 
M is the number of states in Sp, R is the time cost of the function next, 
and E is the time cost of checking syntactic equality of two states. 


Proof: For every state G € Sp, we make exactly two calls to next (line 
3). This yields a factor M2R. The number of iterations is bounded by 
M, since at least one new state must be added in each iteration. Adding 
elements to a list can be done in time proportional to the length of the 
list. The length of Trans is bounded by 2M, and the length of NewStates 
by M. Hence the list operations in lines 4 and 5 can be carried out in 
time O(M). The list NewTrans has length at most 2M, hence computing 
targets(NewTrans) can be done in time O(M) resulting also in a list of 
length at most 27. Removing duplicates from a list of length | requires at 
most 1? comparisons. Hence remdup(targets(NewTrans)) can be computed 
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in time O(EM?). Finally, with a similar argument, removing elements 
from remdup(targets(NewTrans)) that occur in States can also be done in 
time O(EM7?). Summing up, we obtain an overall complexity of O(M2R + 
M(M+M + EM? + EM?)) € O(RM + EM). 


We now relate M, R and E to the length of the initial specification F. 
The numeric interpretation of F, specifically the value K, will be used in 
the complexity analysis. The next lemma relates kK and len(F). 


Lemma 6 For the initial specification F, log(K) € O(len(F)). 


Proof: We first make the following observation. A closed, integral expres- 
sion which maximises value with respect to length is of the form X” for some 
n EN. Since len(X”") = 2n—1 and val(X") = 2” it follows that for all closed, 
integral function expressions C, val(C) < 2!"(), It now follows that 


K= \d| ae |m| fe |n| < plen(Ca) es glen(Cm) 24 glen(Cn) <3: glen(F) 


and hence that log( A’) € O(len(F)). 


Next we show that the length of function expressions in Sp can be 
bounded in terms of the length of F. 


Lemma 7 For all function expressions G € Sz, len(G) € O(len(F)?). 


Proof: We first show that for all closed, integral function expressions C, 
len(C) € O(log(|val(C)|)?). (8) 


A closed, integral C € FExpr which maximises len(C) with respect to |val(C)| 
is of the form C = —(1+X+X?+4+...+X") for some n € N. One can easily 
prove (by induction on n) that len(—(1 +X +X? 4+...4X")) =n? 4+n+2. 
Since n < log(|val(C)|), it follows that len(C) € O(log(|val(C)|)?). 

By the definition of the next-function, any G in Sp is of the form G = 
(D+ (Cm xX s))/Cn, hence len(G) < len(D) + len(F). By (8) and Lemma 4, 
len(D) € O(log(K)?), and hence by Lemma 6, len(D) € O(len(F)?). It 
follows that len(G) € O(len(F)?) + O(len(F)) = O(len(F)?). 


In order to express R in terms of len(F’), we analyse the length of 
expressions produced by the function nextp. As a first step, we show that 
derivatives of a closed, integral expression C are linearly bounded by len(C). 
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Lemma 8 /f C is a closed, integral function expression in polynomial nor- 
mal form and a € 2, then len(Cqa) < 3- len(C). 


Proof: We first show that for all k > 1: 
len(X*) = 6k —5 (9) 


The proof is by straightforward induction on k. For k = 1, len(Xa) = 
len(1) =1=6-1-—5. Now let k > 1. 


xe = (x?) x x), = (RR & x) 4 (0 1). 
Hence 
len(X*®) = len(X*~1) + 6 =1y 6(k — 1) —-5 +6 =6k—5 


To prove the lemma, we first consider the case where val(C) > 0. Since C is 
in polynomial normal form, C is a sum of terms of the form X*, k € N. For 
all k > 1 and a € 2, we have from (9): 


len(X*) = 6k —5 < 6k —3 =3- len(x*). 


From this observation it follows easily that len(C,) < 3- len(C). 
In case val(C) < 0, then C = —E for some E, and C, = —(E, + L(Ela])), 
hence 


len(Ca) = 3+ len(E,) <343- len(E) = 3- (1+ len(E)) = 3- len(C) 


where the inequality follows from the previous case. 


Now we can give a bound on the length of expressions that must be 
reduced to polynomial normal form. 


Lemma 9 For all G € Sp and a € 2, len(nextp(G,a)) € O(len(F)’). 


Proof: Let G = (D+(C x s))/(1+ (X x E)). Writing out the definition of 
nextp(G, a), we find that 


len(nextp(G,a)) <  len(Da) + len(Cq) + len(E)) +7 
(Lemmas) < 3(len(D) + len(C)) + len(E) +7 
€ O(len(G)) 
(Lemma 7) € O(len(F)?). 
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We can now state the time complexity of our synthesis algorithm ex- 
pressed in terms of the length of the initial specification. 


Theorem 9 The synthesis algorithm which constructs from initial specifi- 
cation F the Mealy machine (Sp,mg) runs in time 20(er®)”) , 


Proof: We first express R in terms of len(F). Let G be in Sp. The time 
complexity of computing next(G,a) is dominated by the time complexity 
of computing pnf(nextp(G,a)). Computing the polynomial normal form of 
an integral expression E is in the worst case exponential in len(E) due to 
the duplication of subexpressions when applying the distributive law. E.g. 
Px (Q+R)=PxQ+PxR. All other manipulations in the computation of 
pnf(E) are polynomial in the length of the expression they are applied to. 
By Lemma 9, len(nextp(G,a)) € O(len(F)”), hence 


RE 20(len(F)”) (10) 
Checking whether two expressions in Sp are syntactically equal is linear in 
the length of the two expressions. Hence by Lemma 7, 

E € O(len(F)?) (11) 
Lemma 4 can be used to show that the number of derivatives of [F] is at 
most K = |d| + |m| + |n|, and in many cases an even tighter upper bound 
can be given. A detailed proof of this result can be found in Theorem 3.3.10 
of [6]. Consequently,  < Kk and by Lemma 6, 

M € o(ien®)) (12) 
Combining Proposition 3 with equations (10), (11) and (12) we find that 
the overall complexity of the synthesis algorithm is 


O(RM Eg EM?*) = O(20Cen(F)?) ; glen(F) at len(F)? : Diente) Ec gO(len(F)?) | 


Remark 1 The complexity result stated in Theorem 3.5.10 of [6] differs 
from Theorem 9 above due to the following. Firstly, the synthesis algorithm 
in [6] is based on computing actual syntactic derivativatives rather than the 
function next which only applies to rational function expressions. In par- 
ticular, the algorithm described in [6] also allows the partial construction 
of Mealy machines from non-rational function expressions. Secondly, the 
syntax of function expressions here differs slightly from the 2-adic expres- 
sions in [6] where terms of the form X" are atomic grammar entities. The 
implementation in [7] is based on the algorithm presented in [6]. 
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6 Discussion and Related Work 


Brzozowski [3] showed how to construct a deterministic finite automaton for 
a rational expression by computing its finitely many derivatives, herewith 
lifting the well-known fact that rational languages have a finite number 
of (left) quotients, to the symbolic level of expressions. Since then, vari- 
ous applications and generalisations have been studied. In [1], Antimirov 
introduced the notion of partial derivative and used it to construct non- 
deterministic finite automata. In [14, 16], we reformulated Brzozowski’s 
original approach in coalgebraic terms and generalised it to formal power 
series over arbitrary semirings, providing at the same time a generalisa- 
tion of Antimirov’s results. A similar generalisation to formal power series 
and rational expressions with multiplicities was found, independently, by 
Lombardy and Sakarovitch [12]. 

We have shown how to construct (in exponential quadratic time) a finite 
Mealy machine realisation of rational 2-adic functions by a symbolic com- 
putation of derivatives. The same principles can be used to perform Mealy 
synthesis of functions specified in mod-2 arithmetic (cf. [6, Ch. 3]). We ex- 
pect that the method can also be extended to include bitstream functions 
specified in the Boolean bitstream algebra, respectively Kleene bitstream 
algebra, described in [17] alongside the 2-adic and mod-2 bitstream alge- 
bras. It would be interesting to combine the operators of these bitstream 
algebras into “mixed specifications”. Although a Mealy machine of “mixed 
expressions” can be defined (using the stream differential equations), and 
various “mixed identities” are proved in [17], it is not clear whether there 
exists an equivalence on mixed expressions with finite index, which is ef- 
fectively decidable. Such an equivalence is necessary to generalise symbolic 
synthesis to mixed specifications. 

Closely related to the work presented here is the work in [19] where coal- 
gebras are synthesised from a kind of generalised recursive process specifica- 
tions. The method is similar to ours in the sense that it relies on a symbolic 
computation of generated subcoalgebras. This construction is parametric 
in the functor T which defines the coalgebra type as well as the syntax 
of the specifications, and so covers many different automaton-like system 
types (including Mealy machines) in a uniform framework. Instantiating 
the results of [19] to Mealy machines, we point out the following differences 
with our work. The syntax of the specification language in [19] is in close 
correspondence with the semantic structure and specifying rational 2-adic 
functions would be inconvenient to say the least. On the other hand, the 
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close connection between syntax and semantics ensures that any behaviour 
of a finite Mealy machine can be specified in their language and vice versa, 
all specifications have finite realisations. As we have seen in Example 8, 
rational function expressions are not expressively complete with respect to 
finite Mealy machines. However, function expressions can specify infinite- 
state behaviours such as f(a) = 0 x o which are not expressible in the 
language of [19]. Finally, we mention that the synthesis algorithm in [19] 
does not necessarily produce a minimal Mealy machine as is the case with 
our synthesis algorithm. 

The idea of generating behaviour from syntax is possible at a very 
general level. Such interplay between algebra (syntax) and coalgebra (be- 
haviour) can often be captured by so-called bialgebras for a distributive 
law (cf. [20, 2]). For example, the stream differential equations for the 2- 
adic operators define a distributive law of the 2-adic signature over stream 
behaviour. Other examples of distributive laws are rules in structural op- 
erational semantics (cf. [2]), and also regular expressions and deterministic 
automata form an example of this much more abstract setup (cf. [11]). 
Such a bialgebraic picture also exists for function expressions and Mealy 
machines, although in a slightly less direct way than in the examples just 
mentioned. This result can be found in [9]. 
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Appendix 
Proof of Proposition 1. Let a € 2 anda € 2”. 
Ola = [0](a:c)(0) = 0(0)=0 = Ola 
OJa(o) = [O](a:o)’ = [0 = (0 = [0a](e) 
ijla] -= [alteze)(O) = 100)=1 = “ile 
tao) = [Maco = = =) = Palo) 
X][e] = [X](a:0)(0) = xX(0)=0 = X{a 
XJa(o) = [X](a:o)! — ee, eat = [X ](c) 
s|[a =. ‘([sl(@ia))\(0) =~ (@are)(0) =a .=- ‘sla 
sla(o) = ([sl(a:0)’ = (aio'=o = [sal(0) 
Note that for any b € 2 and o € 2”: [v(b)] = [6]. 
Minus: 
-F][a] == ([-F](a:0))(0) = (-[F](a:¢))(0) = ([FI(a:2))(0) = [FI lal. 
= Fla] = —Fla] 
-FJa(o) = ([-F](a:0))! = (-[F](a:0))' 
= ~(FI(a:0) + (ele: 


Sum: 


— 

=) 

as 
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= > 

32 

<- @ ae aa 

Oo sogee 

Ss: Se = 

“foe 2 

oe +< Oo —~ 
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Ss ee Se 

Bee ee ey 

Re OS a Tak igs AO 
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Soe Se ee 
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Product: 


([F](@:2) x [G](a:2))(0) 


([F x G](a:o))(0) 


[F x G][a] 


([F](a:0) x [G](a:0))’ 


Inverse: 
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Proof of Lemma 3. Let f(co) = ise for integers d,m and n with n 
odd, and let a € 2. First, using the stream differential equations of Defini- 
tion 1 and the identities of commutative rings, we find by a straightforward 
calculation that the initial value and stream derivative of a bitstream quo- 
tient o/7 (with 7(0) = 1) are given by: 


(a/7)(0) =o(0) and (a/r)' = (0' = [o(0)] x 7')/7 (13) 


We now use (13) to compute fa(o) for a € 2. The steps marked with (tf) 
use commutativity of x and +. 


d+(m x eo) (t) (He x my 


n n 


fa(o) = ( 


(d+ ((a:0) x m))' = ([d(0) @ (aA m(0))] xv‘) 


d' + ((ao x m) + ([a] x m’)) + [d(0) A (a A m(0))|— 
([d(0) 6 (aA m(0))] x n’) 


(d! — ({a(0)] x nt) + (mx 0) es 
Ai n 
7 ((d' +m! + [d(0) A m(0)]) — ([d(0) & m(0)] x ne’) 
se ifa=1 


n 
The initial values are computed in a similar manner. Hence we obtain the 
following equations for the 6(a)-value in (5): 


6(0) = d'—[d(0)] x n’ and 
6(1) = d’ +m’ + [d(0) A m(0)] — [d(0) 6 m(0)] x n’. 


The rest of the proof is now straightforward using (2) (p. 104). If d is even 
then 6(0) = d' = $d, and if d is odd, we get: 


| i 1 
= = 1 1) = 
6(0) =d—n 5 (4 ) 5 (n ) 5 (4 n) 
When d+ ™ is even with d and m both odd, then 
ee 1 1 1 
6(1)=d +m +tl=5(d 1)4 5 (m 1)+1=5(d+m). 


If d+ ™m is odd with d odd, and m even, then 


1 1 1 1 
6(1) =d' +m! —n! = <(d 1) +5m 5 (n 1) = 5(d+m—n). 


The remaining cases are proved similarly, details are left to the reader. 


Proof of Lemma 4. It is a consequence of Lemma 3 that the derivatives 
of f have the given format (6), since f is itself of the form required in 
Lemma 3, and hence so are all derivatives of f. We prove by induction 
on the length of w € 2* that the numeric value 6(w) is in the given range. 
The base case (w = ¢€) is clear. To prove the inductive step, we use the 
numeric interpretation of derivatives of rational 2-adic functions given in 
Lemma 3. To ease notation, let 1 = min{d,—n+1,-n+m+1} and u = 
max{d,m-—1,0}. Note that | < 0 < u. Assume as induction hypothesis (IH) 
that 1 < d(w) < wu. Inequalities obtained from the induction hypothesis 
will be denoted by <r; y. 

Induction step for 6(w0): We first consider the case where 6(w) is even, 
and thus 6(w0) = 56(w). We have the following cases: 


Now if 5(w) is odd, then 6(w0) = $(6(w) —n). To prove the lower bound, 
we have 


1 < -n+1 > 21 < l—-n+1 <ry d(w) -—n+1, 


and since 6(w) —n + 1 is odd, it follows that 21 < d(w) — n and hence | < 
$(6(w) —n). The upper bound follows easily from 6(w)—n < 6(w) <rH u 
which implies 5(6(w) —n) < u, since u > 0. 

Induction step for 6(w1): If 6(w)+m is even, then 6(w1) = $(6(w)+m). 


We first prove the lower bound. We have (since n > 0), 
I< -n+m+1<m and I <ry 4d(w) 


whence 21 < 6(w) +m, andl < $(6(w) +m). For the upper bound, we 
have 
m—-1< wu and 6(w)<ryu 


whence 6(w)+m—1 < 2u, and 6(w)+m-—1 must be odd, since 6(w)+m is 
even. It follows that 6(w) +m < 2u, which in turn implies $(6(w)+m) < uw. 
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If 6(w) + m is odd, then 6(w1) = $(d(w) + m—n). We know that 
1 < —n+m+1 and hence 21 < l—-n+m+1 <ryq 6(w)+m—n+t+1. 
Since 6(w) +m—n-+1 is odd, it follows that 21 < 6(w)+m-—vn, and hence 
L < $(6(w) +m-—n). 

The upper bound is proven in a similar fashion. We have m—1 < u 
which implies d(w) +m—n <7yp utm—n < u+m-—1 < 2u, and it 
follows that 5(6(w) +m—n) < u. 
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